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ABSTRACT 


Least squares methods are used with NIMBUS IV SIRS-B 
clear-column radiances as independent variables and two 
identical-period sources of geopotential heights (at the 
500 mb and 300 mb pressure levels) as dependent variables. 
Regression equations are developed for each of the sets of 
geopotential heights for three latitude bands in the 
Northern Hemisphere. Three-day and four-day data bases were 
used on each regression specification. Each set of equa- 
tions was tested on independent data at a subsequent 
composite period of 12-hours and 24-hours following the 
dependent pooled samples. The regression results of both 
the dependent and independent tests were examined for 


possible operational usefulness. 
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I. INTRODUCTION 


The Satellite Infrared Spectrometer (SIRS-B) carried on 
the NIMBUS IV spacecraft system relays to National 
Aeronautics and Space Administration (NASA), at Goddard 
Space Flight Center, radiance observations taken in seven 
distinct channels of the 15 wm band of carbon dioxide, and 


one at 899.0 cm”? 


in the atmospheric window. In addition, 

there are six channels in the water-vapor rotational band 

also being relayed by NIMBUS IV, but these latter six 

channels will not bear on the subject of this dissertation. 
Each of the eight primary channels (699.3 стт, 


L 714.3 cm”, 


677.8 cm +, 692.3 cm, 699.3 cm“, 706.3 cm 
750.0 Er and 899.3 sm is sensitive to temperature 
variations in the tropospheric-stratospheric range of heights. 
In addition, each channel has a maximum transmissivity in a 
discrete height range above the surface. Smith, Woolf, and 
Jacob (1970) presented a graph depicting dT /dlnp versus 
pressure where T is the transmissivity and the derivative, 
dt /dlnp, becomes a weight factor of the relative contribu- 
tions of the various elevations to the total emitted 

radiance in the channel. It can be seen in Figure 1 and 
Equation (1), after Smith et al (1970), that the different 
channel radiances are most responsive to changes in tempera- 


ture which occur at the level of their peaks in dr /dln p. 


Each channel tends to signal the occurrence of temperature 


T2 





changes at or near its own pressure-level peak, as 


discernible from Figure 1. 


PRESSURE [МВ] 





Fig. 1 Derivative of transmittance with respect 
to the logarithm of pressure for eight 
SIRS-A channels [ after Smith et al 
(1970) ] . 

The radiative transfer integral applied to a clear 
column of air for which the sounding T(p) is Known gives 
the radiance as 

(Pg) 
N(V) = Β[νι,τίρο)] TUAE B NAO] ат; (1) 
1 
where B [vi ro] is the Planck radiance at wave number 


v and  T(p), the temperature at pressure level р. In 


1» 





Equation (1) the subscript "s" indicates a surface value and 
τίνι, P) is the fractional transmittance of the atmosphere 
in the spectral interval centered at vj, from pressure 
level p up to the satellite. ТЕ the sounding is known, 
the radiance in the channel centered at wave number У; 
can be computed since τί Ν., p) is known as a function of 
pressure p Smith (1970) . The inverse problem has become 
more relevant since the advent of the SIRS instrument, 
namely in the use of observations in the eight channels to 
deduce the clear atmosphere structure T = T(p). The usual 
method of solving an integral equation for the function 
B( Vi, T) was found to be unstable, but it was found | Smith 
et al (1970) | that regression methods relating temperatures 
and geopotential heights at standard pressure levels were 
quite well described by the results of SIRS scan-soundings 
Smith, Woolf, and Jacob (1970) found it convenient to 
convert the clear column radiances into equivalent 
"brightness temperatures," OV). The definition of 
TB( Vi) is based upon the particular Planck function for 
temperature Tg( Vi) which equals the clear-column channel 


radiance N( Vi), that is 


alv y. Tp vi) | - NCy 1) (2) 
so that by solving for Tg we have 


Tp (V ;) = CoV; E | avv) + 1 | (3) 
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Here the constants Cj and Co have the values 


Cy 1. 19061 š 107? erg-cm?-sec^i -steradian | 


1.43868 cm-°K 


> 
Smith et al (1970) derived regression equations for contour 
heights at standard pressure levels (see Table 1) utilizing 
the eight simultaneous brightness temperatures T&( Vi) as 
independent variables. As dependent variables they employed 
a large sample of radiosonde-reported geopotential-heights 
Ec the mandatory levels. “Their results for the month of 
September, 1969, using a two-week sample, was stated in the 


nonlinear regression form 
8 
203) = 20) + $55 Gppp[T67) - vo] + 


i-1 


(4) 
8 
узву EA = тө (к) 2 


with an analogous regression equation for T(p;). Their 
regression coefficients for the September, 1969 SIRS-A 
example are shown in Table 1, where the coefficients for 
Т(р.) are also listed. For the data period, 24-29 December 
1970, employed in this study, an analogous set of regression 
equations based upon updated SIRS-climatology would be used. 
It is clear that the form in Equation (4), as employed 
by Smith et al (1970), essentially builds from the tempera- 
ture structure of the remotely-sensed atmospheric columns 


with p. = 1000 mb and reconstructs a height profile by 
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statistically relating Z(pj) through the hydrostatic 
equation to "Z(p,.) at any other level desired: 

The objective of this study was to obtain statistical 
relationships for geopotential heights of the 500 mb and 
300 mb pressure levels by stepwise regression analysis 
between these variables and the eight clear-column satellite 
radiances. Further, it was desired to do this without prior 
computation of the brightness temperatures. The latter 
procedure is an unnecessary step if corrected radiances are 
known, as was the case in this study. If a sixteen- 
predictor equation were developed, Equation (4) would then 
be modified to include the М, апа N,* terms (i = 1,...8) 
with appropriate coefficients. 

With the use of NIMBUS IV SIRS-B corrected radiances, 
N;, and National Environmental Satellite Center (NESC) 
regression-computed contour heights at the same geographic 
locations, diagnostic regression equations for the 500 mb 
and 300 mb geopotential heights were developed. Separate 
regression equations were developed using both three-day 
and four-day pooled data bases. The regression equations 
developed were tested for validity on an independent sample 
in each case. 

The same procedure was applied to Fleet Numerical 
Weather Central (FNWC) analyzed geopotential heights at the 
500 mb and 300 mb levels for the same positions and times. 
In the latter case, the contour heights had in no way been 


‘related to SIRS-B scan-soundings. Thus, any diagnostic and 
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prognostic relationships discovered here would have to be 
based upon the statistical correlation coefficients 


R[2(5,3),N,], i= 1,...,8, as previously suggested by 
Smith et al. 
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II. DATA PROCESSING 


A. THE ORIGINAL NESC DATA 

The original data obtained from the National Environmental 
Satellite Center (NESC) consisted of a computer listing, by 
days, of selected SIRS-B scan spots during the period 
Ber 1970 - Jan. 1971. Included in this listing were the 
geographical sites and times of each observation, amounts 
and heights of effective cloud cover, temperatures and 
heights of the standard levels, and a double set of 
radiances. The latter were the original set of uncorrected 
radiances and the cloud-corrected values for each channel.- 
The uncorrected radiances in the window channel were missing 
since that channel had become inoperative by December, 1970. 
The sample of scan-spots presented by NESC was therefore 
restricted to sites having known surface temperatures so 
that a realistic effective window channel radiance, B(«',T&g) 
for a surface pressure level pg = 1000 mb could be assigned 
for each location. 

It should be pointed out at this stage that the SIRS-B 
radiometer carried aboard the NIMBUS IV differed from the 
SIRS-A radiometer aboard the NIMBUS III in that the SIRS-B 
radiometer had a "side-looking" capability as well as the 
vertical mode. This permitted a wider range of soundings 
in SIRS-B than was possible with SIRS-A. These side-view 


radiances vere normalized to a single atmospheric 


E9 





emitter-depth located over the scan spot so that all the 
data was compatible, apart from variations in cloud cover 
and surface-layer lapse variations. 

Because of the scan-spot area of (225 km) 2, most spots 
viewed were partially cloud-filled with some effective 
fraction, A. This fact complicated the determination of the 
cloud-corrected T(p) profiles since the radiances were 


reduced to 


N(V;) =A | B | — Ро | rv 9.) 
f: | 
1 


+ (1 = A) | B Ka T (i. P.) 


/ 


1 


E [| v TO) | T 
(5) 


Ts 


в | Мрт) | а τι) 


assuming a one-level cloud model with the top at p.. Here 
the part of Equation (5) within the second brace is that 
Bart of N(Y;) which represents the desired uncontaminated 
radiance in channel i which escapes through the cloud-free 
fraction, 1 - А. 

Smith et al (1970) have outlined a method for determining 
Br cloud-contamination correct2on which yields the estimate 
for corrected radiance N.(V ;) in each channel. Since the 
work in this study was dependent upon the accuracy of the 


corrected radiances, a brief explanation of this correction 
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procedure is given, based on the one-level cloud model with 


mp at level ре. Equation (5) can be broken =intomtwo parts) 
Noi) = N( Vj) - ε(νι) (6) 


where the correction C(4 i) is designed to eliminate the 
effects of cloud-reduction of the radiance, and may be 
obtained by dividing Equation (5) by the factor l-A. How- 
ever, this will not be a permissible procedure since the 
Cloud cover, A, may be very close to unity. , Instead, Smith 
et al (1970) suggest an iterative procedure whereby м) 


of Equation (6) may be obtained in the form 
COMEDIA | м.) х [ М. .Ре/ ТР) | | (7) 


The function x [i Po: T (p) | represents the contents of the 


second brace of Equation (5). Successive corrections 
c(v;) are then determined for each channel 4,....,8, by 
making the assumption that the radiances observed in 
channels 2 and 3 (750.0 cm} and 714.0 cm 1) are correct 
during iterations while correcting the remaining set of 


radiances Ny. ---Ng. 


In determining a first guess to a clear column atmosphere 
T(p) the five most opaque channels in Figure l are assumed 
correct at the first iteration, as well as the pre- 
determined window channel. The NESC correction scheme then 
makes use of a six-predictor regression equation for  T(p) 


(analogous to Equation (4) but lacking channel 2 and 3 terms) 


DD 





which then specifies the required first guess T (pj) at all 


Standard levels. With this profile now specified the 


best-fit consistent choices of A, l-A, and Po which comes 


closest to providing the best-fit solutions N(4/i) to 
Equation (5) for channels 2 and 3 are obtained. In this 
test p, is subject to variation over the standard levels 
Š = 200, 250,....850 mb and the best-fit A, l-A, are 
thereby chosen. At each iteration, there is an improved set 
of six corrected radiances, which leads to an improved 
specification T(p;). This then leads to improved newly 
Serrected radiances Nel vj) (i = 4,....,8) in accord with 
Equations (6) and (7). Finally, in a limited number of 
iterations C(a⁄ ;) approaches a convergent value. At this 
final stage, the six-predictor regression equation for 
T(p;) gives estimates of corrected radiances in channels 2 
MaA’. The listed statistics relating p, and A are of 
no particular importance since these quantities relate only 
to an effective cloud cover over an area (225 km) 2, and the 
eweud—depth is considered irrelevant in obtaining муб). 
The important factor to be tested in this paper is the 
concept that geopotential heights as taken from a large 
radiosonde sample are highly correlated with the set of 
clear-column radiances, i.e., to the estimated T(p) 
structure of the atmosphere. Thus the NESC SIRS-B data 
obtained by the December 1970 version of the NESC regression 


equation system provided also the geopotential heights and 
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temperatures at all standard levels in the atmosphere from 
p = 1000 mb to 10 mb. 

In developing the NESC regression equations, Smith et al 
(1970) utilized a data sample of approximately two weeks, 
continually updating the sample period closer to the time 
when the independent calculations were made. The heights 
used in their regression equations were interpolated in 
space and time from National Meteorological Center analyzed 
fields at standard pressure levels to the geographical loca- 
tions and observed times of the SIRS measurements. Only 
portions of analyses and NIMBUS III data close to radiosonde 
observations were used for developing their dependent data 
regression coefficients giving heights from corrected 
radiances. These NESC dependent samples were divided into 
latitude bands (18-40°N, 35-55°N, and 50-80°N). A similar 
technique was used by NESC in making available the data from 
December 1970, a portion of which was used in this thesis. 
This data consisted of corrected radiances at scan-spots 
which were used in conjunction with an eight-predictor 
regression equation for generating temperature and height 
values for the standard levels from the SIRS observations. 
These computed heights and temperatures appeared on the NESC 
data listings for 24-28 December 1970, which were used in 
this study. 

The NESC regression-derived heights of the 500 mb and 
300 mb levels, position, time, and clear-column radiances 
for each sounding from 24-27 December, 1970, were used for 


„һе regression analysis described in Chapter III. 
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В. FNWC DATA 

Height fields for 0000 GMT апа 1200 GMT covering the 
period 24-29 December 1970, were obtained from the data- 
tapes of Fleet Numerical Weather Central (FNWC) at Monterey 
for both the 500 mb and 300 mb levels. The contour heights 
were interpolated in space and time to the same positions 
and times of the satellite data obtained from NESC. The 
spatial values of geopotential height were computed by use 
of a double-quadratic Bessel interpolation scheme applied to 
both the synoptic times preceding and following each satel- 
lite σα ο. Then the heights were interpolated 
linearly to the time of the satellite observation. The 
computer program used to accomplish the interpolation scheme 
is listed in the ix. 

The FNWC data provided an independent set of geopotential 
heights for the 500 mb and 300 mb levels from which a second 
set of regression relations could be derived for comparison 


with the regression analyses performed on the NESC data. 
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III. DEPENDENT DATA REGRESSIONS 


A. DEVELOPMENT OF REGRESSION EQUATIONS USING NESC GENERATED 
VALUES OF 500 MB AND 300 MB GEOPOTENTIAL HEIGHTS. 

The NESC-generated 300 mb and 500 mb pressure levels 
(Zy(5) and Zy(3) respectively) were selected for investi- 
gation with the corresponding 500 mb and 300 mb FNWC 
analyses, which presently have no SIRS input. A multiple 
linear regression analysis technique was used in this study 
corresponding to the 8-predictor scheme of Smith, Woolf, and 
Jacob (1970), but containing no nonlinear terms. The speci- 
fic program used was BIMED O2R, a version of the stepwise 
linear regression technique currently on file at the 
Statistical Library of the W. R. Church Computer Facility 
at the Naval Postgraduate School, after Dixon (1966). The 
result of this technique when applied to simultaneous data 
samples of the form [2;, N41. N5, - -, Ng] is to produce 


regression equations of the form 


* 


8 
Z = C, + у CN, = (cC o pnm s : (8) 
=! i 


Here the subscripts 1,..,8 indicate the channel numbering 
system and are also used for variable numbering. A series 


Of dependent test regressions were performed on both  Zw(5) 
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and Zy(3) using a pooled three-day data base and a pooled 
four-day data base of simultaneous values of SIRS-B correc- 
ted radiances. Due to the nature of the data available, 
each data-day presented a composite of scan-soundings 
centered about 1200 GMT, varying up to six hours on either 
Side of this central time. The data was further divided 
into three overlapping latitude bands. The limits of each 


band are listed below: 


Band 1. 209N - 40?N 
Band 2 35ΟΝ - 55ΟΝ 
Band 3 509N - 709N 


The resulting coefficients (Сс. cues) for the three- 
day and four-day regression equations are listed in Tables 2 
and 3 for the 500 mb and 300 mb dependent tests, 


mespectively. 
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TABLE 2, 


Sample Size 


Sample Size 


* 
ludrcates an insignsficonb ΕἘΌΠΈΕΙΡΠΕ Πες τ 


G Variance DY sits soe sre core 


500 mb NESC Dependent Regression 
Coefficients for Three- and Four-Day 
Data Base, Giving Best-fit for Eq. (8) 


3-Day Sample Base 


Band 1 


171 


5473.77 
6.66 


=W 56 
LOS 
28.06 

209722 


=6 60 
2.120 
205 


4-Day Sample 


Band 1 


түт 


5497.55 
6.53 
-15.98 
20713 
34.33 
πο. 
SS 


24.18 
-3.26 
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Band 2 


__ 190 


5689.80 
SO O 


5857 
41.83 
=10.83 
2667/8 


4.02 


19962 
-1.40 


Base 


Band 2 


21 


27250970 
БЕО 
=, 
45.24 

5 
-64.77 
2.64 


19.54 
= 06 


Вапа 3 
177 


4620.48 
=2 790 


= ος 
62.59 
ιν 
UE 


-4.14 
11.74 


Band 3 
216 
4555.32 
-4.45 
-6.38 
54.55 
-10.93 
-35.75 
-7.94 


1.00 
1.45 


analysis 








TABLE 3. 


Sample Size 


Co 


Size 


* 
Indicates an insignificant contribution to 


Best Fit for Eq. 


(8) 


3-Day Sample Base 


Band 1 
teL 


8976.59 
12:79 


=22_ Ὁ 
4.04 


87.54 
-105.20 
5270 


23.90 
5.67 


Band 2 


190 


8490.63 
4.00 


-6.08 
21.62 


47.07 
-96.10 
4.64 


15% SS 
162 


4-Day Sample Base 


Band 1 


206 


8986.32 
13.26 
-21.76 


* 


95.49 
ΤΙ 
-14.34 
28.03 
1. 25 


Band 2 


227 


8613.84 
6.65 
-13.75 
29:09 
41.53 
93.15 
35113 
16.26 
1.65 


of variance by this predictor. 


28 


300 mb NESC Dependent Regression Coefficients 
for Three-and Four-Day Data Base, 


Giving 


Band 3 
177 


7166.52 
1.57 


-15.54 
54.11 


23 l 
-48.47 
-4.34 


6.14 
1.29 


Band 3 
216 


7102.40 


* 


ο. Е 
45. 98 


30792 
-47.49 
-8.48 
ЭЯ 
2.94 


the analysis 








The coefficients in Tables 2 and 3 are considerably 
different in each latitude band which tends to justify the 
use of separate latitude bands in the analyses. Whether the 
limits of each latitude band are optimum or not is undeter- 
mined, but it was felt that any further reduction in the 
width of each band would reduce the number of soundings used 
in each regression equation and thereby have a detrimental 
effect on a statistical analysis. The differences between 
the corresponding coefficients from the three- and four-day 
sample regressions were generally insignificant; however, 
this conclusion will be further —À by examination of the 
fractional explainted variance accounted for by the appro- 
priate regression equations when applied to the data- 
samples. 

In both Tables 2 and 3, the coefficients corresponding 
to channels 3, 4, and 5 appear to be weighted more heavily 
than any other triad of coefficients, which appears to be 
consistent with the response levels of Figure l. 

A useful statistic for the fractional unexplained 
variance is 


2 


1 2 AS k — 1 std. error 
_ u Ὦ - 1 std. dev. (9) 


where k is the number of predictors and n the sample 


size, 
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Thus the fractional explained variance is easily 
obtained from Equation (9) as 


ne ш td 2 | 
R2 = 1 — std. error 
| n l | sta, des. (10) 


Values for R^ and of standard errors for each of the 
dependent tests corresponding to Tables 2 and 3 are listed 


in Table 4. 


TABLE 4. ΠΒ and Standard Errors (gpm) for NESC 
Dependent Regression Equations for 
500 mb and.300 mb Three- and Four-day 





samples. 
500 mb 
3-Day Sample Band 1 Band 2 Band 3 
R? . 8050 .8222 `` |. .9050 
Std. error (m) 61.72 77.13 69.68 
4-Day Sample 
R2 . 8039 . 8154 .9051 
Std. error (m) 60.55 ` 76.78 , 69.79 
300 mb 
3-Day Sample Band 1 Band 2 Band 3 
R2 .8790 .8693 , .9351 
Std. error (m) τι 5] 91.80. 75.36 
u 
4-Day Sample 
R? „8764 „8622 „9349 
Std. error(m) 70222 91.42 75.50 


> 
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As indicated in Table 4, Band 3 appears to produce higher 
values for RZ, and therefore better results, than the other 
two bands in the NESC dependent regression analysis. The 
standard errors for the 300 mb tests were generally higher 
than the 500 mb tests, but one might expect this to happen 
since the standard deviations of Z2,(3) are greater than 
those of .Z2,(5). In general, the 300 mb tests gave higher 
values of κ; upon specification than the corresponding 


500 mb tests. 


B. DEVELOPMENT OF REGRESSION EQUATIONS USING 
FNWC-INTERPOLATED VALUES OF 500 MB AND 300 MB 


GEOPOTENTIAL HEIGHTS. 


Procedures similar to those used in Part A of this 
ES er were used on the FNWC 500 mb and 300 mb geopotential 
heights [Zp(5) and Zp(3) respectively ] , corresponding 
to the identical SIRS scan-sounding positions and times. 
The same latitude band limits as specified in III.A. were 
used for both the three- and four-day regression analyses. 
The resulting coefficients for the FNWC regression equations 
are listed in Tables 5 and 6 for the 500 mb and 300 mb 


dependent tests, respectively. 


Sl 








TABLE 5. 500 mb FNWC Dependent Regression Coefficients 
for Three- and Four-day Data Base 


3-DAY SAMPLE BASE 


WP E de τ; ai 

ample Size 
on 5617.09 5583.51 4615.09 
Cy 7.58 -2.25 -2.50 
Co -19.82 -2.58 -9.43 
C3 23.50 38.35 57.19 
En 25.55 -6.14 -16.06 
Ce -53.38 -67.11 -32.17 
C6 -13.71 5.85 -9.10 
C7 23.14 19.19 Ίο 
ορ -1.35 -2.37 0.80 


4-DAY SAMPLE BASE 


Band 1 Band 2 Band 3 
Sample Size 206 227 216 

Co 5658.75 5762. 27 4595.81 
Cy 7.70 0.36 -2.24 
e -20.43 -10.14 -9.16 
C4 24.21 46.38 55.7 2 
CA 28786 -13.41 -15.03 
5 -59.44 -64.52 -30.50 
C6 -12.29 4.74 -9.58 
Cy 27.24 29312 8:71 
Са -5.01 -2.15 3226 
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TABLE ©. 


Sample Size 


“Indicates an insignificant contribution to 


Size 


3-DAY SAMPLE BASE 


Band 1 
ШУД 


9156.82 
14.03 
-33.45 
28.34 
50799 


-72.67. 


292322 
22.93 
8.82 


Вапа 2 


190 


8451.37 
2:59 9 
=5 2 
34.19 
2923 
266507 
5295 
15.05 
0.78 


4-DAY SAMPLE BASE 


Band 1 


206 


9226.42 
15201 


36.75 
31:79 
52:29 

-77 . 84 

-27.24 
26.96 


3.44 


variance by this predictor. 
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Band 2 
227 


8702.09 


7.09 
220798 
49.50 
16.47 
=83:39 
4.64 
13972 
0.41 


300 mb FNWC Dependent Regression Coefficients 
for Three- and Four-day Data Base 


Band 3 
177 


ἜΘΗ 


κ 


-18.36 
71.49 
26956 

295.292 
πο 6 

9218 


* 


Band 3 
216 


7209825 


* 


=16G 37 
66.94 


* 


οσο 
πω ο 


5.07 
2. 85 


the analysis of 








О “να 


The FNWC regression coefficients listed in Tables 5 and 
6 lend support to the observations from the NESC regression 
coefficients, in that the coefficients from the FNWC 
dependent test vary sufficiently between latitude bands to 
justify using separate bands. Furthermore, the differences 
between the three-day samples and four-day samples are again 
small; and the coefficients corresponding to Channels 3, 4, 
and 5 appear to be weighted more heavily than the others. 

Table 7 lists the values for R2 and standard errors 
for each of the FNWC dependent tests listed in Tables 5 and 
6. As in the NESC tests, the values for R^ were higher in 
Band 3, standard errors were generally larger for the 300 mb 
pests than for the 500 mb tests, and the values for RZ? 
were higher for the 300 mb tests than for the 500 mb tests. 
A comparison of the three- and four-day sample regressions 
for both the 300 mb and 500 mb cases revealed no significant 
differences in corresponding values of R2, when subjected to 


to a t-test. 
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TABLE 7. RZ and Standard Errors (gpm) for FNWC 
Dependent Regression Equations for 


500 mb and 300 mb Three- and Four-day 
Dependent Samples. 


500 mb 
—3-Day Sample Bande] Band 2 Band 3 
R? „7910 „7954 „9014 
Std. error (m) 62.53 84.55 72.48 
4-Day Sample 
R2 er „7899 „8989 
Std. error (m) 63.79 83.69 73.21 
300 mb 
3-Day Sample Band 1 Band 2 Band 3 
R2 .8575 .8361 9235 
Std. error (m) 84.38 105.99 84.32 
4-Day Sample 
R2 „8459 „8328 79222 
Std. error(m) 85%77 104.26 85.07 
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C. COMPARISONS BETWEEN NESC AND FNWC DEPENDENT REGRESSION 
ANALYSES 

If one compares the results from Tables 2, 3, 5, and 6, 
channels 3, 4, and 5 appear to be weighted more heavily than 
the other channels for both the NESC and FNWC dependent 
tests. Table 8 combines the results of Tables 4 and 7 to 
facilitate comparison of NESC and FNWC dependent regression 
test results. Here, only comparative values of R2 by the 
various dependent stratifications are considered. τε 
standard errors are required, reference should be made to 
Tables 4 and 7. 

Band 3 gave the highest values of R^, averaging 8.726 
migher over all dependent regression systems than Bands 1 
and 2, varying from 5.6% to 12.0% higher for individual 
comparisons. The 300 mb values of R2 averaged 4.65% 
higher overall than the 500 mb values of R2. with the NESC 
300 mb values of R2? averaging 5.01% higher than the NESC 
500 mb values, and the FNWC 300 mb values of R2 averaging 
4.38% higher than the FNWC 500 mb values. 

Although quite close in most instances, in all cases the 
NESC values of R were slightly higher man the correspond- 
ing values of R from the FNWC tests, averaging 2.01% 
higher for the NESC tests. A closer Fook reveals an average 
of 1.69% smaller values, or "shrinkage," for the FNWC 500 mb 
tests as compared to the corresponding NESC tests. An 
average shrinkage of 2.32% occurred for the 300 mb FNWC 


‘tests as compared to the NESC tests? tne που щы 
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shrinkages of R^ for each band, averaged over both the three- 
and four-day data-samples is shown for the FNWC tests rela- 


tive to the NESC corresponding tests: 


Band 1 Band 2 Band 3 
500 mb 1.9096 2.6196 0.49% 
300 mb 2.60% 3.13% 1 721 


Bands 3 shows the smallest shrinkage of Re 


Values for FNWC as 
compared to the NESC dependent tests. That the NESC values 
of R? are higher than the FNWC values is not too surprising 
since the NESC data was originally the result of essentially 
an eight-predictor equation, the precise form of which was 
not known in this study. 

A comparison of three- and four-day overall regression 
systems showed that the three-day values of R^ were .46X 
higher than the four-day values. Since this is such a small 
difference, it is likely to be caused only by the additional 
data samples used in converting the particular three-day 
sample to a four-day sample, and the difference was not 
considered significant. 

Thus far, it is seen from the dependent-data tests on 
NESC and FNWC samples that the results are generally good, 
considering the R2 values averaged .8579, with averages of 
.8679 for NESC and .8478 for FNWC. From the paragraphs 
preceding this one, it appears that the dependent tests 


produced better results in Band 3 than the other bands, 
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better results at 300 mb than at 500 mb and better results 
for NESC than for FNWC specification. 

It should be pointed out that the NESC and FNWC height 
data were highly correlated. The 500 mb values of R? ranged 
from .9267 to .9699 and the standard errors varied from 35.26 
meters to 43.15 meters. The values of R2 for the 300 mb case 
ranged from .9244 to .9707 and the standard errors varied 


from 51.46 meters to 59.64 meters. 
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IV. INDEPENDENT SAMPLE TESTS 


A. VERIFICATION PROCEDURE FOR THE INDEPENDENT DATA TESTS 

Fach of the dependent regression equations, stratified by 
NESC and FNWC, by latitude bands, and for the different sizes 
of data-bases in Chapter III, was tested using independent 
data samples. The equations developed from the three-day 
samples were tested with independent data from a fourth day, 
for times centered at 12 hours and 24 hours later. The 
regression equations er from the four-day samples were 
likewise tested for times centered 12 and 24 hours later 
using the equations developed from data collected over 24-27 
December, 1970. The corresponding four-day independent tests 
utilized data from the 28th (28/00 GMT and 28/12 GMT, i.e., 
12 and 24 hours later). In addition, the 12- and 24-hour 
tests provided homogeneous statistics and were pooled to 
provide a larger number of independent data samples for each 
band and height. 

The results of the independent tests generate equations 


of the general form 


21(5,3) = Бо + 012, (5,3) (11) 
апа Zp(5,3) = Бо ЗЕ (522) (12) 


where le) and Zp (5,3) represent the height-estimators 


from the regression equation (Equation 8), D and Fo are 


O 


regression constants, and D, and Г) the regression 


40 





coefficients.  Zp(5,3) are the observed independent-test 
heights determined by interpolation from the data tapes of 
FNWC for the 500 mb and 300 mb levels. All independent 
ZN(5,3) were the result of application of the NESC update 
regression to the independent set of scan-spot corrected 
radiances. Regression analyses for these tests were similar 
to" the -regression techniques applied in Chapter III, in that 
the newly computed height Z is now being compared by 
regression methods to the known height at each scan-sounding 
location. The independent test results are presented in a 
slightly different manner in Tables 9, 10, 11, and 12. The 
mean and standard deviation for the NESC and FNWC verifica- 
tions, and the values for R and standard error resulting 
from Zy(5), Zy(3), Zp(5), and Zp(3) are listed in Tables 9, 


10, 11, and 12 for the three- and four-day independent tests. 
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сли 


Вапа 1 


Вапа 2 


Вапа 3 


Вапа 1 


Вапа 2 


Band 3 


TABLE 9. 


Verification Resulting From the Use of 


Equation (8) Upon 500 mb NESC Independent 


Data. 


Sample 


Size 


12 hir 23 


24 hr 35 
12+24 hr 58 


12 "τ 31 


24 hr 37 
12+24 hr 68 


12 hr 32 
24 hr 39 
12+24 hr 71 


Sample 

size 
12 ΠΤ 28 
24 hr 57 


12+24 hr 85 


12 hr 33 
24 hr 50 
12+24 hr 83 


2 hr 26 
24 hr 47 
12+24 hr 73 


3-DAY TEST 


Mean 


5298526 


5738.60 
5744.41 


5464.45 


5477.38 
5471.48 


525327 
52931 
5289.21 


4-DAY TEST 


Mean 


97500027 
5739158 
5746.36 


5513.99 
5432.87 
5465.11 


5290.45 
5286.16 
5207709 
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Std. 


реу; 


99.04 
121.72 


112.59 


109776 


154.80 
132215 


215-551 
227 2 
220.49 


Std. 


Dev. 
ar piii tta 


90795 
121.50 
Те 


140.14 
164.26 
159237 


202,12 
172.66 
182.30 





7985 


Леб] 
271552 


.5428 


. 1092 
.6934 


+2781 
. 9015 
. 8859 





NE 
‚8308 
8151 


USE 
. 1499 
«1515 


. 7685 
. 82901 
< 7999 


Std. 


Error 


45.52 


57207 
52.64 


69% 30 


76.08 
73.62 


76.54 
72.24 
75.01 


Sta. 


Error 


46.00 
50.55 
48.67 


76.77 
83209 
80.00 


99726 
12238 
82.94 





TABLE 10. 


Band 1 12 
24 
12+24 


Mand 2 12 
24 
12+24 


Band 3 12 


24 
12+24 


Band 1 12 
24 
12+24 


Band 2 12 
24 
12+24 


Band 3 12 
24 
12124 


Verification Resulting from the Use o E 
Equation (8) Upon 300 mb NESC Independent 
Data. 


hr 
hr 
hr 


hr 
hr 
hr 


IE 


hr 
hr 


hr 
hr 
hr 


hr 
Ar 
hr 


hr 
hr 
hr 


Samp 


25 
39 
58 


31 
3 
68 


32 


39 
71 


3-DAY TEST 


le 


Size Mean 


9427.13 
9398.46 


9409.82 - 


8941.29 
8921.95 
8963.41 


8682.44 


8695.46 
8689.59 


4-DAY TEST 


Samp 


28 
57 
85 


33 
50 
83 


26 
47 
73 


le 


Size Mean 


9418.45 
9409.62 
9412.46 


9012.87 
8923.04 
8958.68 


8096.30 
8696 398 
8696.62 
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Std. 


Dev. 


‘166.40 


190790 
180.62 


157.58 
202.08 
122296 


293.24 


303.44 
296.83 


Sed, 


Dev. 


160.38 
203.85 
189270 


206.88 
234.25 
226.81 


266.23 
242.90 


249.60 





. 8980 
. 8526 
. 8658 


71172 
„8027 
‚7704 


‚9084 


. 9334 
29207 





. 8348 
. 8718 
. 8615 


. 7994 
„8137 
«8148 


‚8301 
. 8827 
. 8603 


Sta. 
Error 





54.39 
74.40 
66.76 


85.23 
91.04 
88.34 


90823 


7 9233 
84.19 


Std. 
Error 


66.42 
73.64 
71.08 


94.15 
1027155 
gar 22 


111.299 
84.13 
93.94 





TABLE ll. Verification Resulting From the Use of 


Equation (8) Upon 500 FNWC Independent 


Data o 

3-DAY TEST 
Sample Std. Std. 
Size Mean Dev. R2 Error 
Band 1 La NY 23 5736.63 91.62 „ОР 46.21 
24 hr 35 5724733 130.80 „ТОЗ 71.46 
12+24 hr 58 5729.20 12621] . 7130 62.75 
Band 2 J 2 nr 3] 5468.17 123209 ‚ 6163 77.56 
DIE 37 5473750 159.20 «τ. 82.31 
12+24 hr 68 5471.07 142.86 ‚6950 79.49 
Band 3 ]J2 hr 32 5275.86 208.00 8440 83252 
24 hr 39 5291.72 229758 8863 78.45 
12+24 hr zu 52843257 2918.72 .8627 81.62 

4-DAY TEST 
Sample Std. 2 Std. 
Size Mean Dev. R Error 
Band 1 12 hr 28 57 50287 80.03 . 8329 33.34 
24 Nr 57 5728.56 120204 .8054 53.43 
12+24 hr 85 57352890 108752 .8104 47.53 
Band 2 12 hr 33 5522753 142,47 .6730 83.94 
24 hr 50 542771 171.16 .7143 92.43 
12+24 hr 83 5465.49 166.85 πο 89.31 
Band 3 1; Dr 26 5280.39 198.52 .6612 117.93 
24 Nr 47 5274783 180.50 „ТӨЗ 85.95 
12+24 hr 73 52716.81 185.76 97.26 
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TABLE 12. Verification Resulting From the 
Equation (8) Upon 300 mb FNWC 
Data. 
3-DAY TEST 

Sample Std; 

Size Mean Dev. 
Band 1 Ίος 23 9396.04 160.42 
24 hr 35 9375. 17 195; 55 
12+24 hr 58 οσον ο 161-16 
Band 2 To hr 31 8934.96 187.09 
24 hr 37 8975.52 217555 
12+24 hr 68 8957.03 208.76 
Band 3 12 ας B2 8665.11 282.86 
24 hr oo 8694.04 308.89 
12+24 hr 71 8681.00 295.70 

4-DAY TEST 

Sample Std. 

Size Mean Dev. 
Band 1 I2shr 28 9410.01 144.12 
24 hr 57 9392.98 200.90 
12+24 hr 85 30239219 183.45 
Band 2 τ he ЗЭ 9005.81 217227 
24 hr 50 8913.67 244.64 
12+24 hr 83 8950.30 237. 16 
Band 3 12 hr 26 8078.84 267.39 
24 hr 47 8679.63 257.69 
12+24 hr T3 8679.34 259.33 
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Use of 
Independent 
2 Sta. 
R Error 
. 8465 64.34 
.7643 96.36 
7688 84.10 
. 7281 99.23 
‚’903 107203 
. 407 3 99.04 
8915 94.70 
„9157 90.88 
.9022 93.213 
2 Std. 
R Error 
. 7820 68.57 
8310 83.34 
. 8198 78.34 
.1295 114 ΕΠ 
7718 118.07 
7641 115.91 
. 7858 126.31 
. 8367 104.31 
.8188 111.18 


B. SUMMARY OF INDEPENDENT TEST RESULTS 


1. Results Internal to the Independent Test 


Table 13 lists the combined 12 and 24 hour 
independent-test values of R2 abstracted from Tables 9, 

10, 11, and 12 to facilitate comparisons of test-methods 
contained in those tables. Again our emphasis is on R4 as 
a measure of the efficiency of the verification since this 
statistic has the effect of normalizing the variable standard 
@eviations in Tables 9, 10, 11, апа 12. 

Band 3, as in the dependent tests, gave the highest 
values of RÊ, averaging 6.99% higher than the other bands. 
The 300 mb values of RZ, again, as in the dependent cases, 
were higher than the 500 mb values, averaging 5.83% higher. 
The NESC values of RZ, as compared to the FNWC values of R2, 
wc higher, as before, by an average of 3.60% over all ot 
the independent tests. However, the amount of shrinkage 
between the NESC and FNWC within the independent tests was a 
minimum in Band 2, as compared to the minimum values of 
shrinkage in Band 3 for dependent tests. This phenomenom is 
due to the fact that the independent test of ps fell off 
considerably in Band 2, whereas the corresponding shrinkage 
was somewhat less with FNWC independent tests in Band 2. The 
three-day estimator tests produced average values of RA 
that were .76% higher than the corresponding four-day values 
for R^. This, again, as in the dependent tests, is a small 
difference and is considered to be due to the small sample- 


size, or to the particular data-sample added to obtain the 
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four-day sample, and therefore of little statistical 


significance. 


2. Comparison of Independent- With Dependent-Data 


Specification. 


From the independent tests, the same basic 
conclusions as for the dependent tests are to be drawn since 
the results are quite similar. The independent-test values 
of R2 were generally good with an average of .8004, an 
average of .8184 for R2 . from the independent NESC tests 
and an average of .7824 for the FNWC tests. The highest 
values of R2 were in Band 3, highest at 300 mb, and 
highest for the NESC tests. 

In proceeding from the dependent tests to the 
independent tests, some loss of accuracy is to be expected. 
Values of R2 again will be used for comparisons. Table 
14 is derived from Tables 8 and 13 to present an easier 
comparison of the corresponding values of R^ from the 
dependent tests and the independent tests, with the value 
of shrinkage indicated also. In most cases the shrinkage 
of the 82 values from the FNWC dependent tests to the 
independent tests was slightly larger than the shrinkage 
for the similar NESC tests. The FNWC shrinkage averaged 
1.59% more than the NESC shrinkage. 

The first group of numbers in Table 14 are averaged 
over all bands and over both three- and four-day coeffi- 


cient estimators. In this group, the values of R2 for 
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the NESC 300 mb tests show the lowest relative shrinkage of 
the four shown in the group. This might have been expected 
since the 300 mb NESC test results from both the dependent 

and independent tests showed better results than the others. 

The second group of numbers in Table 14 is averaged 
over the 500 mb and 300 mb level independent tests, and 
over three- and four-day coefficient estimators 2. This 
second group shows relatively small shrinkages in Band l. 
This could possibly be due to the fact that in low latitudes, 
the radiance correction procedure is somewhat more consis- 
tent than in mid-latitudes, and may be affecting the amount 
of shrinkage in a favorable manner. 

The third group of numbers in Table 14 is averaged 
over all three Trace bands, and shows the relative 
shrinkages by use of the Equation (8) based upon the three- 
and four-day dependent coefficient matrices. The amount of 
shrinkage of the three-day independent tests over the four- 
day tests, or vice-versa, is not very obvious, even though 
one might expect that the four-day tests should be more 
representative due to a larger input of data. From the 
findings of this study, it has consistently been determined 
that tests on independent data using the four-day data-base 
estimators are on the average slightly less effective in 
RZ than when the corresponding three-day data base was 


used. 
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TABLE 14. Summary of Dependent to Independent 
Test Shrinkage. 
Test Being Dependent Independent Average 
Compared Value of R2 Value of R2 Shrinkage 
Zy (5) . 8428 . 7878 5.50% 
Zp(5) „8259 „7547 7.12% 
Zn (3) .8929 . 8489 4.40% 
Zp (3) . 8697 „8101 5.96% 
avg. 5.75% 
Zy Band 1 .8411 „8320 „91% 
Zp Band 1 .8183 „7829 3.54% 
Zy Band 2 ‚8423 ¿7572 8.51% 
Zp Band 2 . 8136 „7359 7.77% 
ZN Band 3 NOT 1. . 8657 5.4475 
Zp Band 3 491.16 ‚8283 8.33% 
avg. 5:73% 
Zn (5) 3-day .8441 27/862 5.59% 
Zy (5) 4-day .8415 ‚7873 5.42% 
Zn (3) 3-day . 8945 . 8523 4. 22% 
ZN(3) 4-day -8912 . 8455 4.5776 
Zp(5) 3-day „8293 . 7569 7 . 24% 
Zp(5) 4-day 8225 .7524 7.01% 
2®(3) 3-аау . 8724 8193 5.31% 
Zp(3) 4-day . 8670 . 8009 6.61% 
| avg. 5.75% 





V. CONCLUSIONS 


It has been shown in this study that regression equations 
developed from relatively small data bases of heights and 
corrected SIRS-B radiances can be used to compute heights on 
independent samples with relatively good results. It is felt 
that further testing should be conducted at other levels 
using comparisons of other samples to test the results 
obtained here. If the resources are available, larger data 
bases should be tested in a similar manner. Direct applica- 
tion of similar procedures could be of possible use in the 
update of FNWC analyses for standard levels in areas of 
sparse data, or even aboard ships with weather offices to 
provide real-time height calculations over an operational 
area. In the latter case the advantage of a small data base 
lies in the relatively small amount of core storage necessary 
for the type of computer calculations possible in shipboard 
operations. 

Although the NESC height data provided slightly better 
regression results than the FNWC height data, this was to be 
expected since the NESC heights were computed using an eight- 
predictor regression equation, and therefore were internally 
mere consistent. The results obtained by using FNWC data 
for dependent and independent regression analyses, parallel- 


ing the procedures applied to the NESC data, lends support 


o 





to the overall concept of developing regression equations to 
compute geopotential heights using SIRS-B corrected radiances. 
Improvements of the results should be expected if 
longer-term data samples are used. Improvements should also 
be expected if the width of the latitude bands were smaller, 
and if the radiance data were subjected to gross error checks 
before entering into regression analyses. Also, a fact to 
be considered is the tendency of the SIRS-B channel response 
to degrade and require updated calibration after being in 
orbit for a period of time. If degradation errors were 
linear in time, the regression-climatology should provide 
update corrections, but the normal case is that degradation 
occurs non-uniformly in the different channels. The best 
scan-soundings are ndr ds obtained soon after the satellite 
is launched into its orbit about the earth, before appreci- 


able system-degradation can occur. 
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